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ABSTRACT 



In this paper we consider the process of ahgnment of a spinning black hole and 
a surrounding misaligned accretion disc. We use a simplified set of equations, that 
describe the evolution of the system in the case where the propagation of warping 
disturbances in the accretion disc occurs diffusively, a situation likely to be common 
in the thin discs in Active Galactic Nuclei (AGN). We also allow the direction of the 
hole spin to move under the action of the disc torques. In such a way, the evolution of 
the hole-disc system is computed self-consistently. We consider a number of different 
situations and we explore the relevant parameter range, by varying the location of the 
warp radius i?w and the propagation speed of the warp. We find that the dissipation 
associated with the twisting of the disc results in a large increase in the accretion 
rate through the disc, so that AGN accreting from a misaligned disc are likely to 
be significantly more luminous than those accreting from a flat disc. We compute 
explicitly the time-scales for the warping of the disc and for the alignment process and 
compare our results with earlier estimates based on simplified steady-state solutions. 
We also confirm earlier predictions that, under appropriate circumstances, accretion 
can proceed in a counter-aligned fashion, so that the accreted material will spin-down 
the hole, rather than spinning it up. Our results have implication in a number of 
different observational features of AGN such as the orientation and shape of jets, the 
shape of X-ray iron lines, and the possibility of obscuration and absorption of X-ray 
by the outer disc as well as the general issue of the spin history of black holes during 
their growth. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - galax- 
ies: nuclei 



1 INTRODUCTION 

There are strong observational and theoretical grounds for 
believing that accretion discs around black holes may be 
twisted or warped, in the sense that the plane of the Kep- 
lerian orbits varies continuously with radius and with time. 
From the observational point of view, warped discs have 
been observed in the nuclei of active galaxies (AGN) at rel- 
atively large distances from the central supermassive black 
hole (~ 1 pc), especially th ough water maser obse rvation 
(for exan iple: in NGC4258. iHerrnstein et alj ITool or m 
Circinus, iGreenhiU et "ail 12003) 7 A warped inner accretion 
disc might also be the explanation of the apparent lack 
of correlation between the direction of the radio jets em- 
anating from AGN and the plane of th e host galaxy disc 
jKinnev et al.ll2000l: ISchmitt et 111120021) . Modelling of ac- 
cretion events onto the cores of galaxi es seems to indicate 
that they occur at random orientations iKendall et al .120031: 



ISaitoh fc Wadall2004l) . In addition, a warp in the disc might 
lead to interesting ph enomena of obscuration of the cen- 
tral engine in an AGN JPhinnevllT989i: iGreenhiU et al...2od3: 
iNavakshinlboOSi) . 

From a theoretical perspective, there are a number of 
possible physical mechanisms that can cause the disc to 
be warped. On the large scale, a warp might be produced 
through the interaction of the disc with a smaller mass com- 
panion orbiting in an inclined orbit. However, closer to the 
black hole, at distances of the order of ~ 10^ — 10^i?s, where 
Rs is the Schwarzschild radius of the black hole, the most 
likely cause of warping is due to relativistic effects. In par- 
ticular, if the black hole is spinning and the spin axis of the 
hole is misaligned with the axis of rotation of the disc (or, 
equivalently, if the hole's angular momentum is misaligned 
with respect to the angular momentum of the accreting mat- 
ter), then the Lense-Thirring precession produces a warp in 
the disc. This process is called the Bardeen-Petterson effect 
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jSardeen fc PettersorJ [ToTBII . Since the rate of the Lense- 
Thirring precession is a strong function of radius oc , 
the general result is that the inner parts of the disc tends 
to align the direction (but not necessarily the sense, see be- 
low) of its specific angular momentum with that of the hole, 
while outside a typical r adius, t he warp radiu s, ^ 10^ -Rs 
jNataraian fc Pringlell 998 : Kin g et al.l2005l) . the disc tends 
to retain its original direction of rotation. 

In this paper, we explore the time-dependent evolu- 
tion of the process of alignment between an accretion disc 
and a spinning black hole. We consider the case where 
the p ropagation of the warp oc curs diffusively. This oc- 
curs llPapaloizou fc PringldllQSl when the dimensionless 
viscosity coefficient in the disc (here called Qi, see below) is 
larger than the disc aspect r atio H/R (which should be the 
case for typical AGN discs, IWiiers fc Pringl^ll999ll . Here 
we are not concerned with the detailed physical mecha- 
nism responsible for the propagation of the warp, the discus- 
sion of which would require larg e 3D numerical simulation s 
(c.f. iNelson fc Papaloizoul lioool : iFragile fc AnninosI |2003) . 
Rather, we explore here the long term evolution of the sys- 
tem by using a simplified ID approach (£ringlc 1992). Addi- 
tionally, in our implementation, we allow the black hole spin 
direction to evolve as a consequence of the mutual torques 
with the hole. In this way, the evolution of the hole spin and 
the warping of the disc, which might occur on similar time- 
scales, are computed simultaneously and self-consistently. 
We explicitly evaluate the above-mentioned time-scales, and 
we also determine the effect that the warping process has on 
the accretion phenomenon. In particular, we find that the 
additional dissipation associated with the twisting of the 
disc leads to a large enhancement in the accretion rate (up 
to a factor ~ 10). 

We also address the issue of whether the disc and 
the hole angular momenta might end up being counter- 
aligned. _nus_20SsiMity a,ppea, red ruled out by the anal- 
ysis of IScheue r fc FeileJ lll99d) (hereafter SF). However, 
iKing et alJ l^20^!^^ have recently challenged this conclusion, 
and, using simple geometrical arguments, find that counter- 
alignment should take place when the magnitude of the hole 
angular momentum Jh, the disc angular momentum Jd and 
their initial misalignment are in the following relation: 



^ < -2cos( 

Jh 



(1) 



where 9 is the initial an gle betwe en t he two angular mo- 
menta. [Kin^^^i] (|200^ argue that lSFl did not find counter- 
alignment because they implicitly assume that Jd was in- 
finitely large, therefore making inequality Q impossible to 
be satisfied. Here, we re-examine this issue by performing a 
number of simulations with given Jd and Jh and conclude 
that indeed counter-alignment is possible when inequality 
Q is satisfied. When counter-alignment occurs, accretion of 
matter onto the hole causes the hole to spin-down, rather 
than spin-up. The possibility of counter-alignment is there- 
fore particularly rel evant to the descripti on of the spin his- 
tory of black holes dVolonteri et alj)2005ll and to the deter- 
mination of the innermost stable orbit, which in turn deter- 
mines the efficiency of conversion of matter into luminosity. 

The paper is organized as follows. In Section 2 we de- 
scribe the basic set of equations that we use to determine 
the evolution of the accretion disc and of the spinning black 



hole. In Section 3 we examine more thoroughly the steady- 
state equations found bv lSFt discussing their applicability 
and their limitations. In Section 4 we perform a number of 
simulations describing the evolution of "warped spreading 
rings" , where we consider the evolution of thin rings of mat- 
ter, initially inclined with respect to the hole spin direction, 
as they are accreted and twisted by the black hole. In Section 
5 we discuss our results and draw our conclusions 



2 BASIC EQUATIONS 

In this Section we briefiy describe the model that we use 
to follow the evolution of the system. We consider a sys- 
tem comprising a black hole at i? = (note that here R 
is a purely spherical coordinate and not the cylindrical ra- 
dius), with spin angular momentum Jh, surrounded by a 
disc, with surface density E and specific angular momentum 
density L(ii), in Keplerian rotation around the black hole. 
Since we expect the interesting disc dynamics to occur in the 
region of the warp radius i?w 3> -Rs, treating the black hole 
as a point mass and using Newtonian dynamics is adequate 
for our present purposes. The modulus of the specific angu- 
lar momentum density is thus given by L{R) = E V GM R, 
where AI is the mass of the black hole. The unit vector indi- 
cating the local direction of the specific angular momentum 
in the disc is denoted by \{R). We further denote by Jd the 
total angular momentum of the disc: 



Jd = / 2TvRL{R)dR 



(2) 



We follow the formalism of |Pringi3 il992h and introduce 
the two viscosities vi and 1^2, where vi is the standard shear 
viscosity in a flat disc, and V2 is the viscosity associated 
with vertical shear motions and which describes the diffusion 
of warping disturbances through the disc. The equation of 
motion for th e specific angular momentum in the disc is 
JPringlelll993l : 
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Note that the above equation is a generalization of the 
standard accretion disc equation (e.g. iPringle,.1981i') and re- 
duces to it in the case of a flat disc {81/ dR = 0). The viscos- 
ity coefficients vi and 1/2 , of course, are in general a function 
of local disc properites as well as of R. However, in the rest of 
the paper, we will explore the simple c ase where both vi and 
i/2 are constants, independent on R. IPaoaloizou fc Pringlel 
il983h have shown that in the linear approximation cor- 
responding to small warp angles they are related through 
1^2 /I'l ~ l/2ai, where ai is th e standard a parameter as - 
sociated with the viscosity i/i iShakura fc Sunvaevlll973l . 
In this picture, therefore, the stress associated with warp 
propagation is much larger than the viscous stress, since Qi 
is generally much smaller than unity. However, the actual 
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value of V2 in the non-linear case has still to be determined 
(see discussion in Section 2.1. below). 

Equation Q describes the evolution of a warped disc 
subject only to internal torques. In order to study the 
Bardeen-Petterson effect we add a term on the right-hand 
side of equation @, that describes the relative torque be- 
tween the disc and the hole due to Lense-Thirring precession: 



dL 

dt 



(4) 



where the precession rate fJp is given by (e.g. 
iKumar fc Prineia.lQSa flp = ujp/R^, with: 



2GJh 



(5) 



where c is the speed of light. This is valid for Rs- 

In most of our simulations, we also allow the black hole 
spin to evolve as a consequence of the torques exerted by 
the disc. The equation of motion for the hole spin follows 
directly from Newton's third law and we have: 



dJh 

dt 



(6) 



where the integral is taken over the disc surface. Note that, 
since the torque on the black hole is perpendicular to the 
hole angular momentum, the modulus of Jh is constant with 
time. The only effects of the relative torques are to make the 
hole spin change direction or precess and affect alignment 
with the angular momentum of the disc. The calculations 
here neglect the changes in the magnitude of the spin of the 
hole caused by the accretion of matter onto it. In general 
we expect the accretion induced change in the hole spin to 
occur on a longer time-scale than the disc induced change 
in spin dire c tion both because of the lev er-arm effect at R^ 
jReej|l978t iNataraian fc Pringld Il99il and because warp 
propagation (1/2) typically occurs faster than mass propaga- 
tion (vi). We discuss this issue more in detail below. 

We solve equati ons llsll and @ numerically using the 
scheme described in IPringld li 19921) . We employ a logarith- 
mically spaced radial grid spanning three to four orders of 
magnitude for most simulations. Equation ® is solved using 
a standard fourth-order Runge-Kutta method. 



2.1 Validity of our approach and theoretical 
uncertainties 

Our basic evolutionary equation [Equation Q] assumes that 
the evolution of the disc warp is diffusive and the response 
of the disc to the warping perturbation is completely de- 
termined, in our approach, by the coefficient 1/2- The ba- 
sic motivation for our assumption stems in the analysis of 
[Papaloizou & Pringlc ( 1983), who have shown that the warp 
propagation is indeed diffusive (as long as q > H/R) in the 
linear regime, and assuming that the internal disc motions 
are purely hydrodynamic, and that the viscosity present is 
the standard Navier-Stokes viscosity. As mentioned above, 
FPapaloizo u & Pringlc (1983) also find that 1^2 /vi = l/2ai, 
which is generally much larger than unity. 

In fact, even with these simplifying assumptions, the 
linear, hydrody namic analysis predicts the presence of pre - 
cessional terms dPaoaloizou fc Pringldligil IOgilviell2"003) . 



These are however second order terms and we neglect them 
here. 

However, the behaviour of real discs is surely more 
complex than what is described simply by Equation (|3Jl. 
Turbulence is likely to be controlled by magnetic fields 
iBalbus fc Hawlev^ll998^ . and what kind of viscosity arises 
under these conditions is yet to be fully determined. The 
only attempt at measuri ng 1/2 in the prese nce of MHD tur- 
bulence has been done bv lTorkelsson et alj J200Q'1 . who find, 
within the limitations of a shearing box approximation, that 
also in this case one should expect 1^2 /ui to be relatively 
large. 

Additionally, non-linear effects are going to be impor- 
tant (indeed, as described below, most of our simulations 
deal with non-linear warps). The behaviour of the disc in 
such regimes is much harder to descri be. The only analyti- 
cal t heories of non-lin e ar wa rps are byiMvie (1999, 20q3) 
and lOgilvie fc Dubu'3 ll200lll . These are for purely hydro- 
dynamic discs and they generally assume a polytropic equa- 
tion of state. In addition, such analyses are limited to mildly 
non-linear warps and the resulting non-linear corrections are 
dependent on the det a iled d isc structure (E(i?), H{R), ...). 
In particular, lOgilvid 1^9^ has shown that the first non- 
linear corrections to the diffusion coefficients become impor- 
tant when the amplitude of the warp ijj = R\dl/dR\ > a. 
For these amplitudes, the horizontal shear motions induced 
by the warp become supersonic, inducing s hocks and add i- 
tional dissipation in the disc, not included in lOgilvi3 il99Sr) . 
In the presence of very lar ge warps there is also the possi- 
bility llGammie et al.ll2000ll of additional dissipation caused 
by fluid instabilities. In this case, one possibility is that the 
ratio V2lv\ might be substantially reduced. 

In view of all the above uncertainties, we have opted for 
a straightforward, and somewhat abstract, approach. The 
reason for doing so is that it keeps the physical processes 
involved simple, and does not specialize to any particular 
disc application (such as AGN, X-ray binaries, etc.). This 
enables a basic understanding of the physics involved in the 
alignment (or possibly counter-alignment, see below) pro- 
cess, without adding physically uncertain, and likely unnec- 
essary, complications. 

We have therefore adopted the simple approach of 
Equation (|3J, but allowing the main parameter U2/u\ to lie 
within a large range of values (from very large values, as ex- 
pected from the linear analysis, to smaller values, perhaps 
~ 1, that might arise in strongly non- linear conditions), de- 
termining in this way the dependence of the alignment pro- 
cess on this basic unknown parameter. 

More complex numerical simulations will be needed to 
determine the local response of a thin disc in strongly non- 
linear conditions, but we leave this for future analyses. We 
expect that the behaviour determined here will qualitatively 
describe the disc evolution (using the appropriate values of 
the diffusion coefficient) even in strongly non-linear regimes. 



3 STEADY-STATE SOLUTIONS 

ISFI have found analytical steady-state solutions to equation 
in the presence of the forced precession for small tilt 
angles described in equation Q, subject to various simpli- 
fying assumptions about the disc density distribution and 
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Figure 1. Evolution of the warp in the case where Rw = ujp/v\ = 0.25 (log Hw = —0.6) and U2/V1 = 1 (left panel) and V2/v\ = 10 
(right panel). The lines refer to t = (solid line), t = 0.05 (short-dashed line) and t = 1 (long-dashed line), where the times are expressed 
in units of the viscous time (R? /ui) at _R = 1. The dot-dashed line shows the lSFl solution in this case. At time t, the warped portion of 
the disc extends out to a radius such that t ^ R'^/u2, showing that the warp indeed diffuses out at this time-scale. 




log(t) 
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Figure 2. Evolution of the the x- (left panel) and y-component (right panel) of the torque exerted by the disc on the hole, normalized 
to their asymptotic values in a steady state, in the case where U2/U1 = 10. The x-component, responsible for the alignment of the hole, 
is never very far from the asymptotic value, except in the very early stages of the evolution. On the other hand, the j/-component, 
responsible for the precession of the hole's spin can be significantly different and approaches the steady state value more slowly. 
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the inner boundary. The steady-state warp only depends on 
a single parameter _Rw = ajp/!/2, the warp radius, which 
is the radius at which the warp propagation time-scale 
equals the local forced precession rate. Roughly speaking, 
the disc becomes aligned with the hole spin at R <j; 7?w and 
keeps its original inclination for R ^ _Rw • The ISFl solution 
takes a very compact form in terms of the complex variable 
W = Ix + ily, where Ix and ly are the x and y component 
of the unit vector parallel to the local direction of the disc's 
specific angular momentum, and th e z-a jcis is parallel to the 
black hole spin. In terms of W, the lSFl solution reads: 



-2(1 -j) 



-/?w 



1/2- 



(7) 



where K is a, complex constant that determines the ampli- 
tude of the warp at large radii. This solution is inaccurate 
at very small radii, close to the black hole, and holds only 
under the assumption that the inner disc radius Rin <^ Rv, 
and K <^1. 

Then, ISFl compute the torque exerted by the steady 
state disc on the black hole and are therefore able to compute 
the time-scale over which the hole eventually aligns with the 
outer disc. The surface density distribution and the radial 
profile of the warp are assumed to correspond to that of 
a steadily accreting disc and thus this procedure implicitly 
assumes that the alignment time-scale is much longer than 
the warp propagation time-scale (over which the steady state 
solution is attained). The alignment time-scale is roughly 
given by: 



Jd(-Rw) 1^2 



(8) 



where here Jd(-Rw) is the angular momentum of th e di sc 
within the warp radius. The time-scale computed by ISFl is 
therefore formally valid only if the angular momentum of 
the hole is much larger than the angular momentum of the 
dis c w ithin the warp radius. In the steady-state computed 
bv lSFl it can be easily shown that: 



^align — 3(2 



viM f R 



V2 M V ^■ 



1/2 



(9) 



where a is the spin parameter of the hole, R^ is its 
Schwarzschild radius and M is the steady state accretion 
rate. 

A second l imitat ion of the lSFl solution, as pointed out 
by iKing et alJ (|200^ (see also Introduction), is that their 
steady-state disc formally extends to infinity and there- 
fore has an infinit e total angular momentum. According to 
iKing et alj (|200^, this prevents the possibility of obtain- 
ing a configuration where eventually the disc and the hole 
angular momenta end up being counter-aligned. 

With the abo ve discussion in mind, we now assess the 
validity of the ISFl solution by computing explicitly the evo- 
lution of an initially fiat but tilted disc and following its 
approach to a steady state. In this case, we will keep the 
hole spin direction fixed. The disc extends from _Rin ~ 10~^ 
to -Rout = 10^ . Particular care has been given to the bound- 
ary conditions. In fact, we would like to have a numerical 
scheme that conserves angular momentum as far as possible. 
In order to do this, w e implement th e boundary conditions 
exactly as described in |Pringl3 il992h . in such a way that at 



both the inner and outer boundary no torque is provided to 
the disc in either the azimuthal or the tilt direction (apart 
from the Lense-Thirring torques). Note that the usual inner 
boundary condition adopted in numer ical models of accre - 
tion discs (i.e. E = 0, see, for example, iPringle et al]ll986l) . 
implies a zero viscous torque at the inner edge, thus imply- 
ing an advective loss of angular momentum from the inner 
edge and consequently allowing accretion. In our implemen- 
tation, in order to produce an accreting, but torque-free in- 
ner boundary, we set up a sink of material there. In effect, 
we add a term on the right-and side of Equation ||3J of the 
form: 



9L 

dt 



(10) 



where f{R) is a dimensionless function strongly peaked at 
the inner edge and tsink is an adjustable time-scale, chosen 
to be small enough to approximate the usual E = inner 
boundary condition. 

We feed the disc at a constant mass accretion rate 
AI — 1 between _R = 70 and R — 90, with a flat spatial 
distribution. Time units are determined by setting i^i — 1. 
The initial surface density of the disc is taken to be that 
appropriate to a steady state for a fiat disc in this configu- 
ration: 



Eo 



(i?^/^ - R.^i'')/STTR'^'' Rin <R<R: 



,1/2 



1/2 



(11) 



1/2 



Rll')/3nR 



1/2 



^add < R < R out , 



where 7?add = 80 is the mean radius at which we add matter. 
The disc is initially fiat, and inclined with respect to the 
hole spin by an angle 6. The amplitude of the warp is thus 
given hy K — sin 9. Before running the simulations, we have 
checked that, in the absence of Lense-Thirring precession, 
the disc is indeed in a steady state. 

In this case (where the hole spin is kept fixed), the sys- 
tem only depends on the two parameters v-ijvx and R^. 
Below we discuss separately the two cases where i?„ i?in 
and where _R„ « iiin. Note that we therefore assume that 
7?in > -Rs- 



3.1 Large warp radius 

In this Section we describe the evolution of the system with 
a relatively large warp radius (with respect to the inner disc 
radius) _Rw = ^^jvi ~ 0.25 = 257iin. The initial inclination 
angle between the disc and the hole was here taken to be S = 
0.05. In Fig.0we show the evolution of the local inclination 
between the disc and the hole at different times for the two 
cases where 1^2 /fi ~ 1 (left panel) and 1^2 /fi ~ 10 (right 
panel). The lines refer to t = (solid line), t = 0.05 (short- 
dashed line) and t = 1 (long-dashed line), where the times 
are expressed in units of the viscous time {R? /vi) at R — 1. 
The dot-dashed line shows the lSFl solution in this case. As 
expected, in the case where V2/V1 = 10 the propagation of 
the warp is much faster so that, at any given time, a larger 
portion of the disc is warped in this case with respect to 
the case where 1^2 /v\ = 1- Plotting the same quantity as 
in Fig. Q but where the two simulations are portrayed at 
the same time in units of R? /v2 (the warp diffusion time- 
scale), evaluated at 7? = 1, shows that, in this time units, 
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Figure 3. Relative inclination between the disc and the hole at 
t = 0.5 (solid line), t = 3 (short-dashed line) and at t = 60 (long- 
dashed line) . The steady-state s olut ion in this case is significantly 
different than what predicted bv lSFl (dot-dashed line) . Here i?w = 
0.025 (logiiw = -1.6). 



to infinity. For a disc extending from i?in to i?out the torque 
is given by: 



TsF 



-2-K 



"^"^3 ^ RdR = C(i?i„, i?out)TsF,o, (13) 



where C{Rin, Rout) is a correction factor of order unity, that 
can be easily computed analytically based on the SF solution 
[Eq. CZJ]. 

Fig. |5| shows the evolution of the x- (left panel) and 
y- (right panel) components of the torque, normalized to 
the asymptotic value predicted on the basis of the ISFI so- 
lution, TsF- Initially, the disc is only inclined with respect 
to the hole and does not precess {ly is initially zero), and 
therefore the a;-component of the torque is also vanishing. 
In contrast, the y-component of the torque, dependent on 
the relative inclination and responsible for the precession 
of the hole, is much larger than the asymptotic value. As 
the evolution proceeds, the disc acquires a non-vanishing ly 
and starts to exert an aligning torque, Tx, on the hole. This 
reaches a maximum of ~ I.STsf at t « 0.05, and then slowly 
decreases, approaching the asymptotic value. No te th at the 
value of Tx is never very far from the estimate of lSFl which 
implies that the time-scale given by SF for the alignment of 
the hole is roughly correct. On the other hand, the compo- 
nent Ty of the torque is significantly different than the ISFI 
solution and approaches it much more slowly. This means 
that we expect the hole spin to precess at a significantly 
different and apparently faster rate. 



the warped portion of the disc expands at the same rate in 
the two cases. 

An interesting feature that can be seen from Fig. Q is 
that in the small 1^2 case the inner disc does not p erfectly 
align with the hole, as would be predicted by ISFf, This is 
due to the viscous inward advection of angular momentum 
from the outer disc (neglected in the analysis by SF [cf their 
Equation (6) and the following approximations]. This term 
scales as v\/u2 and indeed the misalignment at small radii 
disappears in the right panel in Fig. where v\/v2 = 0.1. 

An important quantity to be evaluated is the torque 
T = dJh/dt exerted on the hole by the disc. Note that, if we 
take the z-axis to be aligned with the hole's spin, the torque 
has only got the x and y components, which are proportional 
to the y and x components of the disc's angular momentum, 
respectively. Note that it is the x-component of the torque 
which is responsible for the alignment of the hole with the 
disc's angular momentum, while the y-conmonent gives rise 
to a precession of the hole spin. In the ISFl steady state, 
the two components of the torque Tsf.o = {Tsyx, TsFy) are 
equal to each other in magnitude and are given by: 

TsFx = —TsVy^—T^KT.ac\/l-0^V2GM, (12) 

where Eoo is the (almost) constant surface density far from 
the inner edge of the disc. In our case, we have Tioo ~ l/Svr 
(see Eq. IIH . In order to check on which time-scale the 
asymptotic values of the torques expressed in equation l|12|l 
are attained, we have computed the torques on the hole ex- 
erted by our disc in the case where V2/ i'i = 10, for which 
our solution eventually approaches the ISFI solution. Note 
that the values expressed in equation 1121 come from an in- 
tegral over the disc surface, assumed to extend from R = 



3.2 Small warp radius 

Here we describe the evolution of the system with a smaller 
warp radius -Rw = 2.5i?in. The ratio of the two viscosities 
is here taken to be 1^2/^1 = 10. The initial inclination angle 
between the disc and the hole was here taken to be S = 
0.05. Fig.|21describes the evolution of the relative inclination 
between the hole and the disc with time. The solid line refers 
to t = 0.5, the short-dashed line to i = 3, while the long- 
dashed line refers to t = 60 (in units of the viscous time- 
scale at i? = 1) . The disc in this case eventually approaches 
a steady stat e, b ut this is significantly different than what 
predicted bv ISFI (shown in Fig. |3 with a dot-dashed line). 
This is to be expected, since their solution is only valid when 
the warp radius is much greater than the inner disc radius. 
Because the warp radius is not much larger than the inner 
disc radius, the inner disc is flat, but remains misaligned 
with the angular momentum of the hole. 



4 WARPED SPREADING RINGS 

From the discussion of previous section, it appears that an 
initially flat disc approaches the steady-state warped con- 
flguration rather slowly, on a time scale ti/j = R /i^2- The 
surface density of the disc, on the other hand, evolves on 
the viscous time scale t^i = R^/vi- The ratio of the two 
time sc ales is t^^/t^, = v\lv2- Based on the linear calcula- 
tions bv lPapaloizou fc Pringlel il983t) . this ratio is small and 
the warped configuration can be achieved before the disc is 
accreted. However, especially for large warps, the linear cal- 
culations are expected to break down and the two time scales 
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Figure 4. Evolution of the angle between the direction of l(-R) (the local direction of the disc specific angular momentum) and that 
of the black hole angular momentum Jji, at different times during the evolution, for the case logi?w = 0.3 and = Jjj (Note that 
throughout the paper we only plot the relative inclination between the disc and the hole in regions where the disc surface density is not 
negligible). Left panel: the lines refer to i = (solid line), t = 3 X 10~^ (short-dashed line) and 6 X 10"'^ (long-dashed line). Middle 
panel: the lines refer to t = 0.01 (solid line), t = 0.02 (short-dashed line) and t = 0.03 (long-dashed line). Right panel: the lines refer 
to t = 0.04 (solid line), t = 0.05 (short-dashed line) and t = 0.06 (long-dashed line). At t = the matter is in a ring at log J? = 0. By 
t ft! 0.3 the whole disc has aligned with the hole spin. 



might even beco me co mparable. It might then be question- 
able whether the ISFl solution can be actually attained in 
realistic systems. 

In reality, as we have mentioned, the surface density 
profile of the disc is unlikely to be that of a steady disc. 
In practice what we wish to model is an accretion event 
(or a succession of events) onto the central black hole in a 
galactic nucleus. Such events appear to occur at random ori- 
entations relative to the black hole spin. For this reason, in 
order to examine more carefully the dynamical evolution of 
a warped disc under relatively more realistic configurations, 
we have performed a number of other simulations, where we 
evolve a system initially comprising a thin ring of matter 
orbiting around a black hole whose spin is misaligned with 
the direction of the ring rotation. These simulations have 
also the advantage of having a well defined total disc angu- 
lar momentum. In this case, we also allow the hole spin to 
evolve according to equation Q. In this case, we can use a 
slightly smaller dynamical range, since throughout the sim- 
ulations the ring never spreads outside R = 5. Therefore, 
the grid in this case extends from R = lO^'^ to J? = 5 for 
most simulations. In a couple of cases however (when we 
consider small warp radii, see below section 4.2) the inner 
radius was decreased to Rin = 5 x 10~^, in order for it to 
be significantly smaller than the warp radius. We initialize 
the simulations by setting up a thin ring of matter, with a 
Gaussian surface density profile, centered around _Ro = 1 
and with an initial width equal to AR — 10~^. The disc is 
initially flat and inclined by an angle 9 with respect to the 
hole spin. The modulus of the angular momentum of the 
disc is therefore Jd ~ AldVGMRo (where Md is the total 
disc mass) . The governing equations only depend on the pa- 
rameters V2/u\, Jh/Jd, and on the ratio R„/Ro. We have 
performed a number of different simulations by varying the 
parameters above. In the following we will describe in more 
detail the behaviour of the system for a number of different 
initial configurations. 




Figure 5. Evolution of the angle between the hole spin and the z- 
axis for the case Jd = Jh and R^/Rq = 2. If angular momentum 
is conserved, and if the systems becomes eventually completely 
aligned, we would expect the final inclination to be Ri 0.4. The 
lines refer to Vilv\ = 30, 10, 5 and 1 (dotted, solid, short-dashed 
and long-dashed lines, respectively). 



4.1 Alignment 

In this case, we take 7?w ~ 27?o = 200-Rin (so that the ini- 
tial material is deposited just outside the warp radius) and 
1^2 jvi = 10. We also take the initial relative inclination to be 
Q — 7r/4 and the ratio of angular momenta to be Jd/ Jh ~ 1- 
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Figure 6. The case = Jh, Rvj/Ro = 2 and U2 = lOvi. Evo- 
lution of the tip of the unit vector parallel to the black hole spin 
instantaneous direction in the (x',y') plane, perpendicular to the 
direction of the total angular momentum of the system + Jd • 
The hole aligns after roughly half a turn in this plane, showing 
that the precession time-scale is of the same order of the align- 
ment time-scale. 



Figure 7. The case = J^, Rvi/Ro = 2. Evolution of the mass 
accretion rate on the hole in the cases where 1^2 /vi = 30, 10, 5 and 
1 (dotted, solid, short-dashed and long-dashed line, respectively), 
and in the case of a flat spreading ring (dot-dashed line). The 
twisting of the disc and the mixing of angular momenta between 
different radii increases the accretion rate by as much as a factor 
of 3-4. 



In this case, the condition expressed in inequality Q would 
tell us that we expect the disc and the hole to align. 

The evolution of the relative inclination of the disc and 
the hole is shown in Fig. |3 The plots display the angle be- 
tween the local direction of the disc specific angular momen- 
tum 1 and the direction of the hole angular momentum, as a 
function of radius R, at different times during the evolution 
(see caption for details), expressed in units of the viscous 
time-scale at R = Ro, to = Rq/ui. At t — the disc is fiat 
and confined in a narrow ring, but then, as it spreads to 
smaller radii, it also gets warped and the inner disc tends to 
become aligned with the hole. By the time t ~ 0.3to, most of 
the disc has aligned with the hole. This is consistent with the 
propagation of the warp being roughly a factor 1^2/1^1 = 10 
faster than the spreading of the ring. 

The solid line in Fig.|^shows the evolution of the angle 
between the hole spin and the z-axis (which coincides with 
the initial direction of the hole spin). During the evolution of 
the system, angular momentum is not perfectly conserved, 
because mass (and correspondingly angular momentum) is 
being removed from the inner disc, to allow for accretion. 
However, the angular momentum removed in this way is, 
in general, a small fraction of the total angular momentum, 
since the specific angular momentum in the inner disc is a 
factor ^y Rin/Ro ~ 0.1 smaller than the total. Neglecting 
this contribution and assuming that the total angular mo- 
mentum of the system Jh + Jd is conserved, we can get a 
simple relationship between the final direction of the hole 
spin and the initial relative misalignment. In fact, if the sys- 
tem becomes eventually completely aligned, all components 



of the angular momentum should align with the direction of 
the total angular momentum (which is conserved). We then 
get that the angle <f) between the initial and final direction 
of the hole spin is given by: 



1 + ^ cos 6* 



(14) 



where ^ = Jd/Jh- In the present case, where ^ = 1 and 
9 — 7r/4, we expect (j) ~ 0.4, in agreement with the result 
shown in Fig. 13 

During the alignment process, the effect of the torques 
between the hole and the disc is not only to align the relative 
angular momenta, but, as discussed above, also to produce a 
precession of both the angular momentum of different annuli 
in the disc and of the spin of the black hole. In order to 
illustrate this, we plot in Fig. [H] the evolution of the tip 
of the unit vector parallel to the instantaneous direction of 
the black hole spin in a plane {x',y') perpendicular to the 
(constant) direction of the total angular momentum of the 
system. At the end of the alignment process, the hole spin 
should be parallel to the total angular momentum and the 
tip should therefore tend to become very close to the origin 
of this plane. This is indeed the behaviour shown in Fig. 
El The hole spin starts with a relatively large x component 
(it starts on the left in this plot) and, after roughly half a 
turn in this plane, eventually aligns with the total angular 
momentum. This also shows that the precession time-scale is 
of the same order as t he a lignment time-scale in agreement 
with the conclusion of lSFt 

The solid line in Fig.|7|shows the evolution of the mass 
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Figure 8. Evolution of tlie angle between l(-R) and Jjj in the 
case where initially 9 = 3it/4 and Jd/Jh = li log^w = 0.3 and 
= 10. The system eventually becomes counter-aligned The 
lines refer to times t =0, 0.01 and 0.1 in units of T{^^lv\. 



accretion rate with time. As we can see, the presence of 
the warp results in an enhancement of the mass accretion 
rate by as much as a factor 3, with respect to a non- warped 
spreading ri ng calculation (analogous to the standard case 
discussed by |Pringlelll98lh . here shown with a dot-dashed 
line. This is due to the fact that the term associated with 
v-i, induces some energy dissipation in the disc, which then 
is driven more efficiently to s mall ra dii. This process was 
already noted and described bv lPring lc (1992). 

In order to illustrate how the behaviour of the disc 
changes when different values of 1^2 /ui are used, we have 
also run a number of simulations, identical to the previous 
one, but with V2/V1 ~ 30, v^jvx = 5 and V2/i'i ~ 1. In 
actual discs, the ratio 1^2/1^1 might be even larger than our 
largest value, but we expect that the trends derived here will 
hold also for even more extreme viscosity ratios. 

The results of these other simulations are also shown in 
Fig.|5]andEI where the dotted line refers to 1^2/1^1 = 30, the 
short-dashed line to 1/2/1^1 ~ 5 and the long-dashed line to 
^2/1^1 = 1. 

The results shown in Fig. 13 can be easily understood 
in terms of Eqs. ||5J and @, according to which the align- 
ment timescale is inversely proportional to 1^2 /vi- Indeed, 
when a large 1/2 / vi the disc and the hole rapidly get aligned 
between each other. An interesting feature appearing in the 
high V2 lv\ simulations is that the hole inclination overshoots 
somewhat before settling down at the asymptotic value. 

For what concerns the mass accretion rates, as expected, 
the increase in accretion when a warp is present is larger 
for larger values of 1^2 /vi, consistent with our interpretation 
that it is the enhanced dissipation associated with V2 that 
drives most of the accretion process. However, note that the 
enhanced accretion does not scale linearly with V2lv\. 



4.2 Counter-alignment 

As a second test, we have considered the case where, as be- 
fore, -Rwarp = 2i?o, i'2lv\ — 10 and Jd/ J\i = 1. In this case, 
however, we start from initially almost anti-parallel angu- 
lar momenta, with 6 — 37r/4. In this caseJKing et al. (2005[) 
predict that the disc and hole end-up counter-aligned, since 
in this case Jd/Jh = 1 < — 2cos0 ~ 1.41. The evolution of 
the relative inclination of the disc and the hole is shown in 
Fig. IHl The three lines refer to t = (solid line), t = 0.01 
(short-dashed line) and t = 0.1 (long-dashed line). As can 
be seen the spreading of the ring is accompanied by a fast 
counter-alignment of the system (th e angle approaching tt), 
in agreement with lKing et alJ (|200^. 



4.3 Eventual alignment with initial spins almost 
anti-parallel 

As we have seen in the previous Section, during the spread- 
ing of a ring of matter towards the black hole, when the 
initial angular momenta are almost anti-parallel, the inner 
part of the disc becomes rapidly counter-aligned with the 
spin of the hole. In general, the angular momentum of the 
part of the disc that first interacts with the hole is only a 
small fraction of the total. Therefore if the disc and hole are 
initially close to counter-alignment we would expect that 
this material would always counter-align, regardless of the 
tot al angular mome ntum of the disc. However, as noted 
bv lKing et alJ i2005ft . when the disc angular momentum is 
large, we should expect that eventually the system would 
tend to alignment. This brings the question of how this even- 
tual alignment takes place, when initially the inner disc is 
expected to be counter-aligned. In order to answer to this 
question, we have performed another simulation, with the 
following parameters: iiw = 2i?o, V2lv\ = 10 and Q — 37r/4. 
This is similar to the case discussed above. However, we now 
take Jd/Jh = 5. I n this case, Jd/Jh = 5 > — 2cos0 « 1.41, 
and, according to lKing et al.l (|200^, we should have even- 
tual alignment. 

The evolution of the system is shown in Fig. |5] The 
left panels show the evolution of the relative inclination be- 
tween the disc and the hole, while the right panels show the 
corresponding surface density of the disc. As can be seen 
the initial spreading is indeed accompanied by an counter- 
alignment (top panels) . Note the eflect of the enhanced dis- 
sipation due to the warp in the evolution of the surface den- 
sity. In this case more matter flows in with respect to a non- 
warped spreading ring, resulting in a secondary "bump" at 
small radii in the surface density (compare with Fig. 3 of 
|Pringle|[l993 . which shows the same e ffect, but for a smaller 
warp amplitude, and with Fig. 2 of iPringld [T981I which 
shows the evolution of a non- warped disc). The middle pan- 
els show the evolution at later stages. The disc has essen- 
tially broken into two: an inner counter-aligned disc and an 
outer disc, containing most of the disc angular momentum, 
with respect to which the system composed of black hole 
and inner disc tends progressively to align. The inner disc 
becomes smaller in size as time goes by, due to accretion (en- 
hanced by the warping) and due to the fact that the inner 
disc is essentially disconnected from the outer disc. The bot- 
tom plots show the later evolution of the system. The inner 
disc is progressively accreted, while the outer disc becomes 
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Figure 9. Left: Angle between l(-R) and Jjj at various times during the evolution of a simulation where Jd/ Jh = 5 and where the initial 
misalignment is = 3tt/A. Here logRw = 0.3 and V'zjv\ = 10. Top panel : at t = (solid line), t = 0.01 (short-dashed line), t = 0.02 
(long-dashed line). Middle panel: at t = 0.03 (solid line), t = 0.04 (short-dashed line) and t = 0.05 (long-dashed line). Bottom panel: at 
t = 0.0568 (solid line), t = 0.0572 (short-dashed line) and t = 0.0576 (long-dashed line). Time is in units of Rq/ui, where the ring is 
initially at i?o = 1- Right: Corresponding surface densities of the disc. 



Evolution of misaligned discs and black holes 11 




Figure 10. Evolution of tlie mass accretion rate onto the black 
hole for the case shown in Figure |5] where initially almost anti- 
parallel angular momenta eventually align. The dashed line shows 
the corresponding evolution for a fiat spreading ring. The accre- 
tion rate is greatly enhanced during the accretion of the inner 
disc. 



progressively more aligned with the hole. The long-dashed 
line in the bottom plots shows the final stage of the evolu- 
tion, where the inner disc has been completely accreted and 
we are left with an aligned disc. The outer parts of the disc, 
which still have a small inclination with respect to the hole, 
are going to be aligned more slowly, due to the longer time- 
scales in the oute r disc. We then see that the prediction by 
iKing et alJ i2005ft is still valid in this case. 

The evolution of the mass accretion rate onto the hole 
is displayed in Fig. llOl with a solid line, while the dashed line 
shows the corresponding evolution for a flat disc. In this case, 
the strong warping produced in the initial counter-aligned 
configuration leads to a strong enhancement of the accretion 
rate, by factor ~ 6. The accretion rate then suddenly drops 
when the inner disc has been completely accreted and then 
rises again to values of order unity as the final accretion of 
the outer disc takes place. During the first burst of accretion, 
roughly 20% of the disc mass is accreted. 

Figure 1111 shows the evolution of the tip of the unit 
vector parallel to the hole spin in the {x',y') plane, perpen- 
dicular to the direction of the total angular momentum Jt. 
In this case, the initial angle between Jh and Jt is larger 
than 7r/2, even though at later stages the two vectors even- 
tually align. In Fig. llll we display with a dotted line the early 
evolution of the spin direction, during which it has an an- 
gle larger than 7r/2 with respect to Jt, and with a solid line 
the later evolution, when the relative angle is smaller than 
7r/2. The hole spin makes a full turn in this plane, showing 
that also in this case the precession rate is comparable to 
the alignment rate. In the course of the alignment the black 
hole spin axis traces out a large arc on the sky. 



Figure 11. Evolution of the tip of the unit vector parallel to 
the black hole spin instantaneous direction in the (x',y') plane, 
perpendicular to the direction of the total angular momentum of 
the system Jjj -I- Jdi in the case where the initial angular momen- 
tum of the hole is almost anti-parallel with respect to the disc's 
one and where Jd/Jh = 5 (for details see Figure|5J. In this case, 
the hole initially forms an angle larger than tt/2 with Jt. The 
evolution is shown with a dotted line while the angle between Jjj 
and Jt is larger than tt/2, and with a solid line otherwise. Also 
in this case, the precession time-scale is of the same order of the 
alignment time-scale and the alignment takes place after roughly 
one full precession. Note that the black hole spin axis traces out 
a large arc on the sky. 



4.4 Dependence on warp radius 

One important parameter in the dynamics of the system 
considered here is the relative location of the warp radius 
_Rw with respect to the radius within which most of the disc 
angular momentum resides, here identified with Rq. It is not 
simple to estimate the relative location of these two scales 
for realistic models of AG N discs. If we take a typica l AGN 
model, s uch as the one bylCollin-Souffrin fc Dumon3 (Il990l) 
(see also lKing et alJ2005^ . we find that for a typical choice of 
parameters the radius at which the disc angular momentum 
equals the angular momentum of the hole is a few times the 
warp radius. However, one should keep in mind that these 
estimates are based on steady-state models of AGN discs, 
which might not be applicable on the time-scales considered 
here (and are certainly not appropriate in our spreading ring 
calculations) . 

In order to check the effect of the warp radius on the 
alignment process, we have performed a number of other 
simulations, similar to those presented above, but with dif- 
ferent warp radius. In this cases, we have always taken 
u-z/ui = 10. 

First, we considered the case where the initial inclina- 
tion is ^ = 7r/4 and Jd/Jh = 1, in which case we would ex- 
pect the disc to align with the black hole (see section 4.1). 
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However, differently than what discussed in section 4.1, we 
take here _Rw ~ 0.25-Ro, so that the ring is initially com- 
pletely outside the warp radius. Figure fT^ shows the evolu- 
tion of hole inclination (left) and mass accretion rate (right) 
for this case. The dotted lines show the respective results 
for the analogous case with a large warp radius, described 
above in section [4.11 It can be seen that the time-scale for 
the alignment of the hole is longer in the case where the 
warp radius is smaller, in agreement with equation @. The 
peak mass accretion rate on the hole is smaller in this case, 
with respect to the case with a larger warp radius. 

Secondly, we considered the case where the initial in- 
clination is = 37r/4, Jd/Jh ~ 1 and f?w ~ 0.25-Ro- In 
this case, see section 4.2, we expect the disc and the hole 
to counter-align. Also this case does not present significant 
differences with respect to the case with a larger warp ra- 
dius. The disc and the hole counter-align with respect to 
each other. Again, with a smaller warp radius the alignment 
time-scale is longer, as already noted above. 

Thirdly, we considered the more interesting case where 
initially almost anti-parallel spins end up aligned. This could 
be a typical situation in many AGN. We have seen in the 
previous Section that in this case, initially the low angular 
momentum material that first arrives close to the hole tend 
to counter-align with the hole but then, after the accretion 
of this inner counter-aligned disc, the hole and the outer 
disc eventually align. We have performed a number of sim- 
ulations in order to check the dependence of the time spent 
accreting in a counter-aligned fashion on the choice of the 
warp radius 7?w In all these cases the initial inclination is 
6 = 37r/4 and vijvx — 10, but the disc angular momen- 
tum is chosen so as to dominate on the hole, Jd/Jh = 5. 
We have then run a number of simulations with different 
Rvr/Ro ~ 0.25, 0.1 and 0.05. In the last two cases, we have 
decreased the inner disc radius to Rin = 5 x 10~^-Ro, in or- 
der for it to be significantly smaller than the warp radius. 
Fig. ll3l shows the inclination of the inner disc relative to the 
hole as a function of time for the simulation where R^/Rq = 
0.05 (solid line), 0.1 (dotted line), 0.25 (short-dashed line) 
and 2 (long-dashed line). It can be seen that in all cases 
the systems switches from a counter-aligned to an aligned 
configuration, the transition (corresponding to the accretion 
of the inner disc, see previous Section) being progressively 
more gentle as the warp radius decreases. This is because, 
when _Rw is smaller, the alignment between the hole and the 
outer disc is slower [see Eq. Q], so that when the inner disc 
is accreted, the outer disc and the hole are still not aligned. 
As the warp radius decreases, the time spent accreting in an 
counter- aligned fashion increases and it can be a significant 
fraction of the accretion time-scale, up to f ~ 0.2t,^^ (_Ro), 
where t^^ (Ro) is the viscous time at 7i = _Ro. This behaviour 
can be understood qualitatively in the following way. We 
expect that the transition between the counter- aligned and 
the aligned configuration occurs roughly when the angular 
momentum accreted through 7?w becomes comparable with 
that of the hole [cf. Eq. ||5J]. Since the specific angular mo- 
mentum scales as R^^'^, in the case where R^ is smaller, a 
larger amount of matter needs to be accreted before the ac- 
creted angular momentum becomes comparable to that of 
the hole, so that the transition consequently occurs later. 
Equivalently, the same conclusion can be obtained if one 



considers that the alignment time-scale increases with de- 
creasing 7?w [Eq. J^]. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have considered the evolution of a system 
comprising an accretion disc and a central spinning black 
hole, in the case where the hole spin is misaligned with 
the rotation direction of the disc, and the disc is therefore 
subject to Lense-Thirring precession and becomes twisted 
an d warped. Our approach is similar to the one considered 
bv |Pringi3 il992l) . in that we consider the case where the 
warp propaga t es in t he disc in a diffusive manner. However, 
while IPringld (^Q^) considered the warping of an other- 
wise steady-state disc, here we also consider the case where 
the surface density of the disc evolves as the disc is being 
warped. Furthermore, in our evolution scheme, the direction 
of the spin of the hole is allowed to move under the influence 
of the torques exerted by the disc. 

We have carried out a series of simulations under differ- 
ent conditions. First, we have considered the warping of a 
disc which is initially characterized by a steady-state surface 
density distribution (|Pringlelll99l and where the direction 
of the hole spin is kept fixed and considered various cases, 
where the location of the warp radius i?w and the ratio of 
the two viscosities 1^2 /ui are changed. We find that the disc 
eventually settles down in a steady warped con figuration 
(similar to the one obtained analytically bv ISI<1) . but this 
occurs over the relatively long time-scale associated with 
the warp propagation time-scale « R^ /v2, that can be a 
significant fraction of the accretion time-scale. We find that 
if the warp radius is small and the ratio of the two viscosities 
is of order unity, at small radii the warp structure can be 
significantly different from that predicted bv.SFi. 

We also find that the time-scale on which the black hole 
spin and the disc align can be comparable with the propa- 
gation time for the warp. We therefore also considered the 
probably more realistic situation in which a ring of matter 
is accreted by the hole (as can happen, for example, during 
an accretion event onto a central black hole in a galactic 
nucleus). As the ring spreads, it becomes warped and the 
direction of the hole spin is modified. A general feature of 
these warped spreading ring simulations is that the dissi- 
pation associated with the twisting of the disc leads to an 
enhanced accretion rate onto the hole (by as much as a fac- 
tor 10). 

We have computed the detailed evolution of the direc- 
tion of the spin of the hole during the accretion and align- 
ment process. We have found, in agreement with ISfI that 
the "precession" time-scale (i.e. the time-scale of the change 
of the direction of the spin perpendicular to the total angular 
momentum) is comparable to the alignment time-scale, so 
that during alignment the spin of the hole only "goes round" 
at most once (see Figs. 7 and 11). If the jets emanating from 
AGN point in the instantaneous direction of the hole spin, 
our conclusions have important consequences on the shape 
of the observed "p recessing" jets in these systems. Recently, 
iLister et alj i2003tl have reported observations of a "precess- 
ing" jet most likely emanating from an AGN. They have 
fitted the shape of the observed jet by using a model where 
the nozzle precesses slowly around a fixed axis. However, if 
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Figure 12. Left: Evolution of the inclination of the hole spin in the case where the initial misalignment is 9 = 7r/4, J^i/Jh = 1 and 
i?w = 0.25i?o (solid line) and Rw = 2Ho (dashed line). In the small warp radius case, the evolution of the hole spin is slower [cf. equation 
Right; Mass accretion rate on the hole for the two cases. 



the perturbation of the axis of the nozzle is indeed due to 
changes in the spin orientation of the black hole induced by 
the Bardeen-Petterson effect, then our results indicate that 
the assumption of precession a round a fixed axis might not 
be correct. On the other hand. HTister et all i2003ll estimate 
that the typical length-scale of the jet is « 280 pc, and that 
the advance speed of the lobes is « 0.3c. This would imply 
a "precession" time-scale of « 3000 yrs, which is too small 
to be attributed to the Bardeen-Petterson effect [cf. our Eq. 
J^]. In this particular case, the precession is more likely to 
result from spin-orbit coupling in a binary black hole system. 

We would like to stress that the conclusion that align- 
ment time-scale and "precession" time-scale are of the same 
order for the Bardeen-Petterson effect implies that (i) if 
proper precession around a fixed axis is observed, then it is 
most likely not due to Bardeen-Petterson alignment between 
the disc and the black hole and (ii) it provides an observa- 
tional tool to estimate the alignment time-scale, by looking 
at the "precession" time-scale, for those systems where we 
are indeed looking at the Bardeen-Petterson effect. 

We have also considered in detail the issue of the con- 
ditions under which the disc and the hole can be counter- 
aligne d. We have confirmed the suggestion by iKing et alJ 
l|2005|) that, if the initial misalignment between the disc 
and the hole is larger than 90 degrees and the hole angu- 
lar momentum dominates the system, accretion proceeds in 
a counter-aligned fashion. However, also in the case where 
the disc angular momentum dominates and eventual align- 
ment occurs, a large fraction of the mass is initially accreted 
in a counter-aligned fashion. We have shown (see Fig. I13II 
that the amount of time spent accreting in a counter-aligned 
fashion increases as the warp radius decreases. Indeed, the 
alignment time-scale [cf. Eq. (|5J] scales as R^^^^. In a typical 
AGN, the alignment time-scale c an be obt ained combining 
our equation ||^ with Eq. (22) of iKing eTal. (2005) (to es- 
timate R^/Rs). This gives: 




Figure 13. Inclination between the inner disc and the hole as a 
function of time (in units of the viscous time-scale B? /u\ at R = 
Ro for the cases where the initial spins are almost anti-parallel, 
the angular momentum ratio is Jd/Jh = 5 a^nd Rvi/Ro = 0.05 
(solid line), 0.1 (dotted line), 0.25 (short-dashed line) and 2 (long- 
dashed line). 
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where e is the efRciency of accretion, L is the AGN lumi- 
nosity, Le is the Eddington luminosity, Ms is the black hole 
mass in units of 10* M0, and ai and 0:2 are the a parameters 
corresponding to i^i and 1^2, respectively. We therefore see 
that the alignment time-scale is of the order of the lifetime 
of a typical AGN (« 10^ yrs), and we might expect AGN to 
spend a relatively long time accreting in a counter-aligned 
fashion. This counter-aligned accretion mode can have im- 
portant consequences on observables of AGN, such as the 
shape of X-ray iron lines, since for counter-aligned accretion 
the last stable orbit is at a significantly larger radius than 
for the aligned case. 

However, we have also shown that the disc can suffer an 
abrupt change in orientation at a radius close to the warp 
radius, especially during counter-aligned accretion (see Figs. 
4 and 9). A large fraction of the solid angle as seen from 
the center will be filled by the disc itself. This can lead to 
significant absorpti on of X-rays p roduced in the inner disc 
from the outer disc iPhinnevll989|). The surface density of a 
stead y-state AGN disc at R„, is JcoUin-Souffrin fc Dumontj 
|1990D : 

E ^ 2.8 X lO'a-'/' (^) ( -J—] Ml'' 

loos) las) (16) 

which is high enough to obscure visibility of the central ob- 
ject. This implies that direct observation of counter-aligned 
accretion in observationally biased against, so that most of 
those AGN for which the central regions are visible are most 
likely to be co-rotating. 

Our results are also relevant in order to determine the 
black hole spin history during the earlier phases (which 
appear to occur at redshift 1 < z < 5) when the black 
hole acquires most of its mass. During this phase, each 
accretion event increases the black h ole mass significantly 
(AM/M ~ l-3. IVolonteri et al.l2005t) . In general, accretion 
of counter- aligned material spins the hole down. However, 
since the mass accreted during each event is of the order of 
the hole's mass, if all the accreted material rotates about the 
same axis, we would be in a situation like our Fig. 9 and the 
black hole w ould end up aligned and maximally spinning, as 
suggested bv lVolonteri et al] (120051) . However, if this is not 
the case, we might regard one accretion event as a sequence 
of smaller accretion events where at every time the angular 
momentum of the accreted material is small compared to 
the hole. In such a situation, most of the accretion can oc- 
cur in a counter-aligned way and the hole need not end up 
maximally rotating. Thus t he result would be a spread of 
black hole spins (cf Fig. 8 in lVolonteri et alJl2005^ . We con- 
clude that the outcome for the spin history during this phase 
depends crucially on the details of the merging process. 

We should mention that in our approach we have ne- 
glected the evolution of the magnitude of the spin of the 
hole as a consequence of accretion. This, however, is a small 

fraction of the total angular momentum of the system, since 

1/2 

it scales approximately as M^ccR^I^ , where Mace is the mass 
accreted by the hole and Rin- We should also mention that in 
our simulations we have assumed that -Rin = 0.005 — 0.01-Ro, 
where Ro is the initial radius at which matter in injected. 
In actual AGN discs, the inner disc radius (which is of the 



order of the Schwarzschild radius) might be much smaller 
than this, but this does not affect our conclusions as long as 
-Rw is much larger than Rin- 

A further important uncertainty in what we have de- 
scribed in this paper is the magnitude of the viscosity 1^2 
associated with the propagation of the warp. This deter- 
mines the location of th e warp radius and the time -scale 
of the alignment process. IPaoaloizou fc Pringld |l983") have 
shown that in the linear regime 1^2/1^1 = l/2af, where qi is 
the standard viscosity parameter. However, it is not yet clear 
what would happen in the strongly non-linear case described 
here. Numerical simulations up to now have only addressed 
the case where the disc is relatively t hick and warp propa- 
gation occurs in a wave-like manne r llNelson fc Papaloizoul 
ll99gll200d:lFragile fc Anniiiosl2005|) . One possib ilitv is that 
the flow becomes unstable' lGammie' et alj|2000t) . such that 
perhaps 1/2 ~ I'l- In some of the cases discussed here the disc 
can actually break in the sense that the tilt angle shows a 
sharp jump at some radius (see our Fig. 9, for example). 
In such conditions ( see also the numerical simulations by 
iLarwood et alj[l996fl the diffusive equations adopted here 
might not be adequate. To resolve this issue we clearly need 
detailed numerical simulations of warp propagation in thin 
and viscous discs. 
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